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Objectives 

This  research  was  concerned  with  the  dual  problems  of  identifying  the  global 
optimal  structure  and  molecular  aggregates  as  well  as  determining  the  critical  variables 
defining  the  optimum  structure. 

Final  Report  on  Accomplishments 

Substantial  progress  was  made  in  both  the  global  optimization  and  sensitivity 
analysis  of  molecular  structures.  This  is  briefly  summarized  in  the  following: 

1 .  A  deterministic  global  optimization  approach  was  introduced  for  dipeptides  and 
for  oligopeptides  [1,2].  This  approach  combined  the  atomistic  modeling  force 
field  ECEPP/3  with  the  global  optimization  method  aBB.  The  minimization  of 
the  total  potential  energy  was  based  on  the  dihedral  angles  and  consisted  of  a 
novel  branch  and  bound  scheme  in  which  convex  valid  underestimators  were 
derived.  The  proposed  approach  was  shown  to  converge  to  the  global  optimum 
with  theoretical  guarantee. 

2.  A  branch  and  bound  global  optimization  method,  aBB  for  general  constraint 
NLPs  was  proposed  [3].  The  nonlinear  terms  were  classified  as  (1)  bilinear, 
fractional,  signomial  for  which  special  underestimators  were  derived,  and  (ii) 
general  nonconvex  terms  for  which  novel  convex  underestimators  that  employ  a 
novel  DC  transformation  were  derived.  The  proposed  approach  was  shown  to 
converge  to  the  global  optimum  in  a  finite  number  of  steps,  and  extensive 
computational  experience  with  a  variety  of  applications  demonstrated  its 
potential. 

3.  A  novel  theoretical  approach  was  proposed  for  the  generation  of  valid  convex 
underestimators  for  general  twice-differentiable  problems  [4].  This  approach  was 
based  on  the  generating  interval  Hessian  matrices.  Rigorous  lower  bounds  were 
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proposed  for  obtaining  rigorous  eigenvalue  bounds.  Application  of  these 
bounding  techniques  to  highly  nonconvex  problems  that  arise  in  molecular 
modeling  indicated  that  the  obtained  bounds  are  reasonably  tight  to  achieve 
convergence  to  the  global  solution  within  reasonable  computational  effort. 

4.  A  global  optimization  method  was  described  for  identifying  the  global  minimum 
energy  conformation  as  well  as  lower  and  upper  bounds  on  the  global  minimum 
conformer  of  solvated  peptides  [5].  Potential  energy  contributions  were 
calculated  using  the  ECEPP/3  force  field  model.  In  considering  the  effects  of 
hydration,  two  implicit  free  energy  models  were  compared.  One  method  was 
based  on  the  calculation  of  solvent-accessible  surface  areas,  while  the  other  used 
information  on  the  solvent-accessible  volume  of  hydration  shells.  Detailed 
information  on  the  potential  and  solvation  energy  contributions  was  presented  for 
the  terminally  blocked  single  residue  peptides.  In  addition,  based  on  a  procedure 
that  allows  the  exclusion  of  domains  of  the  (4>,<p)  space,  a  number  of  oligopeptide 
structure  prediction  problems  were  considered,  and  the  role  of  the  solvation  model 
in  defining  global  minimum  conformations  was  addressed. 

5.  A  computationally  efficient  algorithm  to  calculate  the  energy  flux  in  a  large 
multibodied  system  was  developed  [6],  along  with  a  visualization  program 
SHOWFLOW.  This  program  was  interfaced  with  the  existing  molecular 
dynamics  packages,  and  results  were  presented  for  the  energy  flow  in  myoglobin 
upon  photolytic  decarboxylation. 

6.  A  molecular  dynamics  study  in  the  active  site  of  myoglobin  was  presented  [7]. 

The  fluxes  were  calculated  from  the  energy  density  derivatives,  while  the 
photolytic  decarboxylation  was  simulated  by  placing  the  carbon  monoxide  group 
on  an  excited-state  potential  energy  surface  until  the  molecule  had  reached  a 
certain  distance  from  the  hemeiron.  An  ensemble  of  27  runs  were  average  to 
distinguish  concerted  energy  flows  from  random  fluctuations. 
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